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ABSTRACT 


>  f  r>  ••? 


t 


For  the  traction  boundary  value  problem  in  nonlinear  elastostatics 
for  a  body  which  is  convex  in  its  undeformed  reference  state  and  with  the 
assumption  of  sufficiently  small  strains  (but  not  necessarily  small  dis¬ 
placement  gradients),  an  upper  bound  is  obtained  for  the  elastic  strain 


energy  in  terms  of  the  L^integral  norms  of  the  surface  tractions  and  body 
forces  with  the  constant  depending  only  upon  the  ratio  of  the  outer  and 
inner  diameters  and  the  physical  constants  of  the  material. 

This  result  extends  previous  known  results  in  linear  elasticity 
(infinitesimal  displacement  gradients)  and  finite  elasticity  (small  but 
finite  displacement  gradients)  into  the  small  strain  theory  of  nonlinear 
elasticity. 
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SIGNIFICANCE  AND  EXPLANATION 


Suppose  that  an  elastic  body  is  in  a  natural  stress-free  state  and  is 
then  deformed  by  the  action  of  applied  surface  tractions  and  internal  body 
forces.  This  situation  is  described  mathematically  by  a  set  of  highly 
complex  nonlinear  partial  differential  equations. 

In  order  to  simplify  the  equations,  it  is  often  assumed  that,  with  the 
possible  exception  of  a  rigid  body  motion,  each  point  in  the  undeformed 
body  has  been  displaced  only  an  infinitesimal  amount  (linear  elasticity)  or, 
more  generally,  by  a  finite  but  small  amount  (finite  elasticity). 

However,  there  are  many  situations  in  which  points  in  the  body  undergo 
relatively  large  displacements,  all  of  which  can  not  be  included  in  a  single 
rigid  motion  of  the  body,  and  yet  the  strains  produced  are  relatively  small. 
This  is  true  for  example  in  the  case  of  long  thin  bodies. 

Therefore,  it  is  desirable  to  obtain  results  using  only  the  assumption 
that  the  strains  are  small  (but  finite)  while  the  displacements  and  displace¬ 
ment  gradients  are  relatively  large.  In  the  present  research,  an  upper 
bound  for  the  energy  in  the  body  in  terms  of  the  surface  tractions  and  body 
forces  is  obtained  using  only  the  small  strain  assumption  whereas  previous 
similar  results  had  been  obtained  under  the  assumptions  of  linear  or  finite 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC  and  not  with  the  author  of  this  report. 


AN  INTEGRAL  BOUND  ON  THE  STRAIN  ENERGY  FOR  THE  TRACTION 
PROBLEM  IN  NONLINEAR  ELASTICITY  WITH  SHALL  STRAINS 

Joseph  J.  Roseman* 

1.  INTRODUCTION. 

An  elastic  isotropic  and  homogeneous  body  is  assumed  to  possess  a  natural  stress- 
free  state  in  which  it  occupies  a  convex  domain  V C  E^ .  It  is  then  subjected  to  the 
combined  action  of  applied  surface  tractions  and  internal  body  forces  and  arrives  at 
a  new  equilibrium  state  in  a  domain  V.  Using  standard  equations  of  nonlinear  elasto- 
statics  and  working  in  the  context  of  small  strain  theory,  we  shall  derive  an  estimate 
for  the  strain  energy  in  the  deformed  body  in  terms  of  the  L^  integral  norms  of  the 
surface  tractions  and  the  body  forces  with  the  constants  involved  depending  upon  the 
material  properties  and  upon  the  geometry  of  the  domain.  The  estimate  is  similar  to 
one  previously  obtained  by  Bramble  and  Payne  [1]  in  the  context  of  linear  elasticity 
and  by  Breuer  and  Roseman  [2]  in  the  context  of  finite  elasticity. 

The  theory  of  mathematical  elastostatics  (cf.  [3],  [41)  assumes  that  the  deforma¬ 
tion  from  0R  to  V  can  be  described  mathematically  as  a  sm< oth  one  to  one  mapping 

from  the  domain  V  onto  the  domain  V  and  that  there  is  a  relation  between  the 
R 

geometric  strain  associated  with  the  mapping  and  the  stress  in  the  deformed  state. 

The  fact  that  the  deformed  body  is  in  equilibrium  inposes  conditions  on  the  stress 
tensor  which  in  turn  leads  to  a  set  of  second  order  nonlinear  elliptic  partial  differ¬ 
ential  equations  for  the  mapping  vector  function.  These  equations,  together  with  the 
known  conditions  at  the  boundary,  constitute  an  elliptic  boundary  value  problem. 

The  equations  of  classical  linear  elastostatics  are  obtained  by  a  linearization 
of  this  nonlinear  system,  a  linearization  which  is,  in  effect,  equivalent  to  treating 
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the  derivatives  of  the  displacement  vector  (the  displacement  being  the  different 
between  the  position  of  a  point  in  V  and  its  forme: r  position  in  1?^)  as  infinitesi¬ 
mally  small  quantities.  In  theory,  therefore,  the  equations  of  linear  elasticity 
describe  the  physical  model  only  for  infinitesimally  small  motions,  although  experience 
has  shown  that  they  give  good  results  in  many  applicat ions. 

In  the  theory  of  finite  elasticity,  it  is  assumed  that  the  displacement  gradients 
are  sufficiently  small,  depending  upon  the  geometry  of  the  body  and  the  material  con¬ 
stants,  but  not  infinitesimal.  Finite  elasticity  should  therefore  provide  a  more  accu¬ 
rate  mathematical  description  of  the  physical  model  than  linear  elasticity  and  results 
in  finite  elasticity  theory  concerning  existence,  uniqueness  and  continuous  dependence 
on  the  given  data  would,  at  the  least,  be  valid  wherever  linear  elasticity  results  are 
valid. 

However,  as  Fritz  John  (5 1  points  out,  in  practice  it  is  possible  to  have  small 
strains  and  relatively  large  displacement  gradients,  especially  in  thin  rods  or  shells 
(indicating  large  relative  "rotation”  of  different  parts  of  the  body).  Thus,  it  would 
be  especially  meaningful  to  obtain  results  based  on  the  nonlinear  equations  assuming 
only  that  the  strains  are  sufficiently  small  (but  not  infinitesimal),  while  the  dis¬ 
placement  gradients  are  bounded  but  are  permitted  to  be  large  in  comparison  with  the 
strains.  That  is  what  is  done  in  this  paper. 

The  techniques  employed  here  are  a  refinement  of  those  used  in  |2]  and  can  be 
extended  in  a  straightforward  manner  to  non-isotropic  and  non-homogeneous  bodies.  The 
convexity  requirement  on  is  imposed  in  order  to  ensure  the  applicability  of  some 

work  of  F.  John  161,  (71  on  the  relationship  between  rotation  and  strain  and  can  be 
relaxed  somewhat. 

Other  results  concerning  bounds  for  the  strain  energy  in  elastostatics  in  terms  of 
the  given  data  include  (11,  [81,  [91,  (101  in  linear  elasticity,  [21,  [111,  [121  in  finite 
elasticity,  and  [13]  for  the  displacement  boundary  value  problem  in  nonlinear  elasticity 
with  small  strains.*  A  description  of  some  of  the  above  mentioned  results  is  given  in  [141. 

*In  both  1131  and  (141,  the  results  in  [131  are  described  as  being  within  the  framework 
of  finite  elasticity.  However,  the  analysis  in  [13J  actually  requires  only  an  a  priori 
assumption  of  small  strains  within  the  elastic  body. 
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2.  PRELIMINARIES . 


In  the  absence  of  surface  tractions  or  internal  body  forces,  an  isotropic,  homo¬ 
geneous  elastic  body  occupies  a  convex  domain  C  with  inner  and  outer  diameters 
d  and  D  respectively  and  with  a  boundary,  which  is  at  least  smooth  enough  for 

the  application  of  the  divergence  theorem.  Under  the  action  of  a  system  of  known  body 
forces  and  surface  tractions,  is  mapped  onto  a  domain  V  with  boundary  jV  with 

a  point  (x  ,x  ,x  )  in  V  being  taken  to  the  point  (y  ,y  ,y  )  in  V.  After  defor- 
mation,  the  body  is  in  static  equilibrium  and  the  mapping  from  to  V  is  assumed 

to  be  one  to  one  and  sufficiently  smootn. 

The  displacement  vector  is  defined  as 


and,  since  the  mapping  is  one  to  one,  may  be  regarded  as  a  function  of  either 

(Vx2,x3)  or  (yx  <y2 <y3>  • 

We  next  define* 
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ik 
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ik 
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ij  kj  ki 
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(2.2) 

(2.3) 

(2.4) 

(2.5) 


ik 

The  tensors  p^  and  P  are,  respectively,  the  covariant  and  contravanant 
tensors,  gik  and  glk  are  the  covariant  and  contravariant  metric  tensors,  and 
is  the  Kronecker  delta  (identity  tensor). 


^Tensor  notation  is  used  throughout  the  paper  and  all  indices  may  take  on  the  values 
1,  2  and  3.  The  summation  convention  is  followed  so  that  a  repeated  index  in  any  terra 
is  summed  over  all  values  of  the  index.  The  magnitudes  of  a  (real)  vector  v^  and 

=  ✓  v.vT  and  |b|  =  / j  b.kb.k. 


matrix 


bik 


are  defined  here  as 


|V| 
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We  have  the  relations 


and 


i  j  jk 

p  p  =  p .  .P  = 

^  1  }k  ik 


i  j  jk 

9  qjk  =  qijq  =  \k 


(.  .  '  •  > 
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There  are  several  definitions  in  the  literature  for  the  strain  tensor  in  nonlinear 
elasticity.  One  of  the  more  common,  which  we  shall  use  here,  is 


eik =  j  (qik 


*..>  . 

ik 


(2.9) 


From  (2.2)  and  (2.4),  this  is  equivalent  to 


i  nu  ■  :5u. 

>  =  I  +  _Jl 

ik  2  '3x 


Xk 


3u ,  3u . 

.  — 1  i 

3x.  3x.  3x, 


(2.9) 


or,  in  terms  of  y.  coordinates, 
i 

1  3u^  3uk  9u^  3^  3u_.  3u^  3u^ 

eik  2  ^3y  +  3y.  Pmk  3y  3y .  Pmi  3y  3y .  PjiPmk  3y .  3y  ^ 

'k  7m  73  7m  7]  J  7j  7m 


(2.19) 


The  medium  xs  assumed  to  be  hyperelastic,  i.e.  it  possesses  a  positive  definite 
strain  energy  density  function  W  =  W(e)  which  has  the  properties  that  for  every 
G  C  V  whose  undeformed  preimage  is  GR  C  the  strain  energy  in  G,  U (G) ,  is 


given  by 


U(G)  =  lU  W (e  Jdx^dx^dx^ 


(2.11) 


and 


W(e,  =  ueikeik  *  +  0(eh  , 


(2.12) 


where  O(e^)  is  a  smooth  term  of  order  of  magnitude  |e|^  with  a  derivative  of  order 
2 

|e|  for  small  |e|  and  X  and  u  are  positive  material  constants. 

In  addition,  there  exist  material  constants  m  and  M  such  that 


m  |  e  |  ^  *•  W(e)  •  M  |  e  | 2 


(2.13) 


if  |e|  is  sufficiently  small. 
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The  stress  at  any  point  in  V  is  represented  by  the  Cauchy  stress  tensor  : 
which  has  the  property  that  for  any  differential  surface  element  ds  in  V  with  unit 
normal  vector  v^,  the  components  of  force  on  this  element  are  TikVs- 

The  body  in  the  deformed  state  is  in  equilibrium  under  the  action  of  the  body 
forces  per  unit  volume  and  surface  tractions  (force/area)  T^. 

The  tensor  satisfies  the  equations 


3t 

=  -f.. 


y^  e  V  , 


and  the  boundary  condi tions 


(2.14a) 


T..n,  *  T.i  y .  e  30  , 
ik  k  1  1 


(2.14b) 


where  n.  is  the  outward  unit  normal  to  30. 

l 

The  existence  of  the  strain  energy  density  function  W(e)  implies  the  following 

relation  between  T.,  and  e.,: 

lk  lk 


ik  det(p)  pkm  3p .  ’ 

lm 


(cf.  Stoker  (3,  pp.  10-14))  or,  from  (2.4)  and  (2.8), 


Tik  det(p)  ^in^km  3e 


(2.15b) 


Equation  (2.15b)  shows  that  the  tensor  is  symmetric  (because  e  is 

symmetric) . 

For  the  problem  considered  here,  it  is  assumed  a  priori  that,  throughout  0,  e' 
is  sufficiently  small  with  respect  to  the  geometry  of  PR  and  the  material  constants; 


more  precisely  we  assume  that 


!e  I 


with  6  sufficiently  small  depending  upon  the  material. 

However,  the  only  restrictions  that  we  require  on  the  displacement  gradients  ar 


that  the  mapping  from  PR  to  0  be  one-to-one,  that  the  magnitudes  of  the  matrices 
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i  )c 

^ik'  **  '  ^ik  an<^  9  and  their  determinants  be  less  than  or  equal  to  two  in 

absolute  value,  and  that  V  satisfy  some  mild  geometric  requirements,  which  are  stated 
in  the  theorem  below. 


3 .  THE  ENERGY  BOUND. 

We  now  state  the  main  theorem:  Consider  an  elastic,  isotropic  body  which  possesses 

a  stress-free  reference  state  in  which  it  is  homogeneous  and  occupies  a  convex  domain 

PR  C  with  inner  and  outer  diameters  d  and  D  respectively.  Suppose  that  the  body 

is  subjected  to  internal  and  external  forces  and  reaches  a  state  of  static  equilibrium 

in  which  it  occupies  a  domain  V  '  F?  with  internal  bodv  forces  F.  (force/unit  volume) 

1 

throughout  V  and  surface  tractions  T.  {force/unit  area)  on  >P,  the  boundarv  of  .n. 

i 

Assume  that  each  point  (x^,x^,x^)  in  is  mapped  to  a  point  (y^  #yn  ,y ^ )  in 

P  in  a  smooth,  one  to  one  manner,  that  the  magnitudes  of  the  covariant  and  contra- 

variant  Jacobian  and  metric  matrices  and  the  determinants  of  these  matrices  are  bounded 

in  absolute  value  at  every  point  by  two,  and  that  the  strain  matrix  satisfies  (2.16). 

Suppose  that  '*P  has  continuous  curvature  with  maximum  principal  curvaturi  not 

greater  than  32/d  and  that  at  every  point  on  -P  there  exists  a  sphere  of  radius  d/04 

which  is  tangent  to  the  surface  at  the  point  and  whose  interior  lies  entirely  within 

P.  In  addition,  there  exists  a  point  in  P,  which  is  taken  to  lx?  the  origin,  with  th* 

properties  that  the  interior  of  the  sphere  )  :  =  —  is  contained  in  0  and 

y .  n  ./Ay ,  y,  >  -  ~  for  all  (y,  rV^,7^)  on  *V ,  where  n.  is  the  outward  unit  normal 
li  kk—  4  12  3  l 

vector . 

Assume  that  there  exists  a  strain  energy  density  function  for  the  material  which 
satisfies  the  relations  (2.12)  and  (2.13). 

Then,  if  4  (Equation  (2.16))  is  sufficiently  small  with  respect  to  the  material 
constants,  the  total  elastic  energy  of  the  body  U(P)  satisfies  the  inequality 

U<0)  1  BdD(|)q[||F||2  +  i||T||2]  ,  (3.1) 

where 

i)  ||  F  ||  is  the  integral  norm  of  |  F  |  over  0, 

ii)  ||t|[  is  the  L  integral  norm  of  |t|  over  30, 

iii)  q  is  a  sufficiently  large  positive  universal  constant,  and 

iv)  B  is  a  constant  which  depends  only  upon  the  physical  constants  of  the 

material. 
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The  rest  of  Section  3  consists  of  the  proof  of  this  result.  Throughout  ‘ . 
following,  all  constants  denoted  by  are  universal,  while  those  denoted  by  n. 

depend  also  upon  the  material  constants. 

We  begin  by  noting  that,  because  of  (2.12),  (2.13)  and  (2.15),  and  because  of  t:.*: 
fact  that  the  magnitudes  and  determinants  of  P,  p  g,  and  g  "  are  uniformly 
bounded,  that  the  integrals  over  either  PR  or  V  of  the  expressions  W(e), 

:..  ,  and  T .  e .  are  all  of  the  same  order  of  magnitude  in  the  sense  that  any  on* 

IK  IK  IK  IK 

is  bounded  in  terms  of  any  other  with  the  constant  depending  only  upon  the  material 
constants.  We  also  see  that 

M  iV(f)  •  (3*2) 


We  now  call  attention  to  the  fact  that  if  the  deformed  body  were  now  subjected  to 
any  rigid  body  motion  (that  is,  any  combination  of  translations  and  rotations),  the 
magnitude  of  the  strain  tensor  e  at  any  point  as  it  moves  with  the  body  would  remain 
invariant  and  so  would  the  magnitude  of  the  stress  tensor  t,  the  strain  energy  V  (P)  , 
and  the  integral  norms  of  the  body  forces  and  the  surface  tractions  (cf.  [15,  pp. 
127-123]).  Therefore,  since  all  the  quantities  of  interest  in  (3.1)  are  unaffected, 
we  shall  assume  that  a  rigid  motion  has  already  taken  place  and  has  been  included  in 
the  mapping  from  V  to  PR.  The  rigid  motion  has  been  chosen  so  that  the  body  is  in 
what  we  consider  to  be  a  convenient  orientation,  as  will  be  explained  below. 

The  papers  of  Fritz  John  T6] ,  [7] ,  on  the  relationship  between  rotation  and  strain 
imply  that  if  the  quantity  <5  in  (2.16)  is  sufficiently  small,  the  body  may  be  rigidly 
rotated  to  a  position  in  which  an  inequality  of  the  following  type  is  valid: 


dx„  dx„dx,  <  K. 
1  2  3—1 


3  W  eikeikdXldX2dX3  • 


d  V, 


(3.3) 


Such  a  position  is  one  to  which  we  have  referred  as  convenient  in  the  paragraph  above. 
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We  note  that  the  integrand  on  the  right  side  of  (3.3)  is  invariant  under  all 


lations  and  rotations  while  the  integrand  on  the  left  depends  upon  the  orientation 
the  body . 

Because  of  the  bounds  on  p  and  p  ^ ,  we  also  have 

3 


lyT  3^7  dYldy2dy3  -  K2  TT  ///  eikeikdyidy2dy3' 
V  k  Jk  d  V 


(2.  A' 


and,  from  the  remarks  in  the  paragraph  preceding  (3.2), 
3u .  ou .  3 

i  l 


l\!  *7  dyldY2dy3  -  K3  TT  ///  Tikeikdyidy2dy3 

V  k  J  k  d  V 


(3.5) 


From  this  point  on,  we  shall  use  dV  ,  dS  ,  dV  and  dS  to  denote  differential 

x  x  y  y 

volume  and  surface  elements  in  V  and  V. 

R 

From  (2.10)  and  from  the  fact  that  is  a  symmetric  tensor,  it  follows  that 


3u. 


<ij  ■  <1/  'ik  a» 


3u .  3u . 

fffp  T  — 1  —  dV 

prok  ik  3y  3y.  y 
V  m  3 


(3.6) 


3u  3u 


1  r  r  r  n  n 

2  'U  PjiPmkTik  3y.  3y 

V  j  m 


dV  . 

y 


Using  the  divergence  theorem  and  noting  (2.14),  we  obtain 

///  TikeikdVv  =  ///  FiUidVy  +  //  VidSv 

pXKixy  piiy  9piiy 


3u .  3u . 

+  f  I  [  p  t  — 1  — i-  dV 

JLJ  pmk  ik  3y  3y .  y 
V  m  j 


(3.7) 


,  3u  3u 

i-fffp  p  T  — £L  _JL  dv  . 
2  PjiPmk  ik  3y.  3ym  Vy 
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Now,  we  have,  from  (3.2)  and  (3.5),  that 


(3.14) 


///  F.u.dV  +  U  T.u.dS  <  r / / /  F.F.dV  ]  ~  T  l / /  v.v.dV  1 

11  y  Ijj  i  i  y  -  i  i  yJ  ?  4  1 

♦  [//  T  T  dS  ]?[//  v  v  ds  ]'  . 

30  1  1  y  30  1  1  V 
Combination  of  (3.14)  and  (3.10)  yields 

ill  T.,«.,dV  <-  2  [  f //  F.F.dV  v.v.dV  ]  ~ 

lk  ik  y  —  u  p  ii  y  i  i  yJ 

+  2  [  [[  T.T.dS  1  ’[//  v.v.dS  . 

L  i  it  yJ  Lvb  i  i  yJ 


We  now  recall  from  the  statement  of  the  theorem  that  y .  n . /  /  y,  y.  ■  ~  ~r  on  ">P 

ii  kk  —  4 

and  that  the  interior  of  the  sphere  j  :  y.y.  =  —  d^  is  contained  in  0.  Next  we 
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refer  to  the  work  of  Bramble  and  Payne  [1]  who  proved  that  if  the  coefficients  a 

and  6.  are  chosen  such  that 
i 


H  v.  dS  =  //  (v.y.  -  v.y.JdS  =  0 

V  i  y  V  4  3  ii  y 


then,  for  a  sufficiently  large  universal  constant  q^. 


//  v.v.dS  <•  Cd  [j]  1  HI  w(  )dv  , 
in  ii  y  ~  5  v1  1 W 


(3.17a) 


^v1 


///  v.v.dV  «  C6d^)  *  ///  W 


VI(fl)dV 


ik  2  3yk  +  3y.^ 


W(8)  =  ),e.k6ik  +  J  xe.  ,ekk  . 


The  quantities  9ik  and  W(9)  represent,  respectively,  the  linearized  strain  and  the 
corresponding  linearized  energy  density  function.  We  note  that 


WCO)  1  <U  +  -  X)H.k9.k  . 


From  (3.12),  we  see  that 


,  3u,  3u 

„  1  f  i  ki 

9.  .  =  r-  [q —  +  F — )  i 

lk  2  'By,  By. 

k  l 
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and  then,  from  (3.20)  and  (2.10), 


,  3u .  3u .  .  ’-u  ,  -u  u 

Q  _  1_  3  _ i.  £  _ _ k  1_  ^  _ n  _ n 

ik  ik  2  ^mk  3y  .  2  mi  3y  ry  2  ^ji^mk  3v .  'V 

m  3  *  3  *  m 


The  elements  of  p. .  and  -  are  uniformly  bounded  and  the  tensor  - 

id  ;y. 

satisfies  (3.5).  Therefore, 


///  eik9ikdVy  -  C7  ^3  ///  Tikeikd 


dV 


and,  from  (3.19), 


///  W  ( 0 )  dV  <  C  ~  ///  re. 

p  y  “  8  a3  0  llc  1 


.  dV 
k  y 


With  q2  =  q  +  3,  Equations  (3.23)  and  (3.17)  imply  that 


q 

//  v.v.dS  <  C  d§)  2  ///  t  e  dV  , 
ii  y  —  9  '■d;  ik  ik  y 


(3. 


and 


///  v.v.dV  <  C  d2  (~)  2  l[ j  x  e  dV 
i  i  y  —  10  '■d-1  ik  ik  y 


p  i  x  y 
Combining  (3.24)  with  (3.15),  we  obtain 


(3. 


V1 

CD)  2 


d(j)  2  ||F||+2^dC 


D-,  2 
9"  '•a-1 


(3. 


Finally,  the  remarks  in  the  paragraph  preceding  (3.2)  imply  that  there  exists 
constant  such  that 


U  (P)  <  C. ,  HI  t,  e ,  dV  , 
-  11  J  y  ik  ik  y 


<3 


•  34a) 


24b) 


25) 


.26) 


and  then  (3.25)  and  (3.26)  lead  to  the  desired  result,  Equation  (3.1). 


4.  SOME  ADDITIONAL  REMARKS. 


As  stated  earlier  in  Section  1,  the  techniques  used  here  are  also  applicable  to 
non-isotropic  and  non -homogeneous  bodies  and  the  restriction  here  to  a  homogeneous 
isotropic  body  was  only  for  the  purpose  of  simplifying  the  presentation.  The  convexity 
assumption  on  can  also  be  relaxed  somewhat;  the  important  thing  is  that  the  domain 

be  such  that  it  is  possible  to  obtain  an  inequality  of  the  form  (3.3)  by  an  appropriate 
rotation.  (See  16],  (7]  for  a  discussion  of  the  factors  involved). 
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